Numerical study of critical properties and hidden orders in dimerized spin ladders 
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Dimerized antiferromagnetic spin-^ ladders are known to exhibit a quantum critical phase transi- 
tion in the ground state, the existence or absence of which is dependent on the dimerization pattern 
of the ladder. The gapped phases cannot be distinguished by the conventional Landau long-range 
order parameter. However, they possess a non-local (hidden) string order parameter, which is non- 
zero in one phase and vanishes in the other. We use an exact diagonalization technique to calculate 
ground state energies, energy gaps and string order parameters of dimerized two- and three-leg 
Heisenberg ladders, as well as a critical scaling analysis to yield estimates of the critical exponents 
v and p. 

PACS numbers: 75.10.Kt, 75.10.Pq, 75.10.Dg 



I. INTRODUCTION 

Spin ladders have been the focus of significant theo- 
retical interest in the past decade or more. This work 
has been grounded by the experimental discovery of sev- 
eral spin ladder and spin-Pcicrls compounds exhibiting 
various fascinating phenomena, 1 Many of the important 
features of these materials have long been understood 
(For reviews see e.g. Refs.[l|0). However, significant dis- 
coveries may yet be made as crucial questions about the 
underlying nature of these systems remain. 

One of the most peculiar properties of spin ladders is 
that the existence of a gap (i.e. mass) depends on the 
number of legs. The spin excitations in a m-leg spin lad- 
der are gapped if m is even, and the system is gapless 
(quantum critical) when the number of legs m is odd. 
The even-m-leg ladders are interesting examples of spin 
liquids, which are currently under enormous scrutiny. 
The spin liquid has, by definition, a gap which is not 
due to long-range order or apparent symmetry breaking. 
This is hard to understand in terms of the basic pos- 
tulates of the conventional Landau paradigm for phase 
transitions,'? 

The purpose of this study is to explore the question 
of the existence of "hidden orders" associated with the 
massive quantum phases of dimerized Heisenberg ladders 
with antiferromagnetic couplings. It is known that single 
dimerized Heisenberg spin-^ chains are gapped, and that 
when two chains are coupled into a ladder, the system is 
gapped even without dimerization^ Counterintuitively, 
the dimerized two- and three-leg ladders can also be gap- 
less. The critical (gapless) lines of the dimerized two- 
and three-leg ladders were first conjectured in Ref. and 
have been confirmed by subsequent numerical and ana- 
lytical work£~— The characterization of the unconven- 
tional quantum phases in such systems in terms of more 
exotic types of orders may ultimately help to uncover 
novel physics which lies beyond the Landau paradigm. 

We present numerical results on finite size spin-^ sys- 
tems with total number of spins (N) of up to N = 30. 



These results are obtained here using an exact diagonal- 
ization (ED) technique based on the Lanczos algorithm. 13 
The total available computer memory is the most impor- 
tant factor limiting the maximum system sizes accessi- 
ble using this method. This is because the dimension 
of the Hamiltonian matrix grows exponentially fast (as 
2 N for a spin-i state). Although restricted to relatively 
small system sizes, the ED method is considered an "un- 
biased method" useful for calculating almost any physical 
quantity. Our implementation of the Lanczos algorithm 
uses periodic boundary conditions and includes an ex- 
plicit restart step with greater than 5 restarts to mitigate 
loss of orthogonality in the Lanczos basis. 



II. TWO-LEG LADDER 



We begin by presenting the results for a two-legged 
ladder with intrinsic dimerization. A similar program is 
carried out for the three-leg ladder, detailed in Sec. III. 
In line with the previous wor k 12 : 14 which used the bond- 
mean-field theory^ we have considered the two possi- 
ble dimerization patterns of the ladder (staggered and 
columnar) on an equal footing (cf. Fig. [1]). 



(a) 



(b) 



FIG. 1: Dimerized two-leg ladder. Bold/thin/dashed lines 
represent the stronger /weaker chain coupling J(l ± 5) and 
rung coupling J±, respectively. Dimerization patterns: (a) - 
staggered; (b)- columnar. 
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The Hamiltonian of the dimerized spin ladder with ra~ 
legs is given by: 

N m 

H = ^^J tt (n)S Q (n).S Q (n + l) 

71=1 a— 1 

TV m— 1 

+ J±EE S Q (n)-S a+1 (n), (1) 

n— 1 a— 1 

where the dimerization occurs along the chains (a = 
1 . . . to) only, with the rung coupling J± taken as con- 
stant. All the spin exchange couplings are antiferromag- 
nctic. The dimerization patterns are then defined as: 

J a (n) = J[l + (-l) n+a 6] (staggered) (2) 
J a (n) = J[l + {-l) n S} (columnar) (3) 

with periodic boundary conditions along the chains and 
open boundary conditions in the rung directions. 
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FIG. 3: (Color online) The ground state energy gap for the 
staggered and columnar two-leg ladder calculated by exact 
diagonalization for N = 28. 



A. Ground State Energies and Gaps 
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FIG. 2: (Color online) The ground state energies for the stag- 
gered and columnar two-leg ladder calculated by exact diag- 
onalization for N = 28. 

We find that the ladders with the columnar dimeriza- 
tion order always have the lower ground state energy, and 
that they are always gapped. The results for the ground 
state energies of the ladders with N = 24 are shown in 
Figure [5] are very close to earlier predictions for these 
systems from the mean-field theory^ 

Also in agreement with the earlier mean-field 
predictions^ our ED calculations confirm that the exis- 
tence or absence of the quantum critical phase transition 
between gapped phases depends on the dimerization pat- 
tern. The columnar arrangement of the strong and weak 
couplings on the legs is never critical, whereas the stag- 
gered configuration possesses a quantum critical point. 



The characteristic results for the singlet-triplet gaps are 
shown in Fig. [3] 

To take into account the finite size effects in the gaps 
obtained by ED and use this data to extract the critical 
indices near the transition point, we have applied a finite 
size scaling analysis^ If we denote the critical coupling 
where the gap vanishes as J± c (the phase diagram of 
the dimerized two-leg ladder is shown in Fig. @| and the 
reduced coupling as g = J± J^ ±C , the gap Ai,(g) in the 
system of the linear dimension L can be expressed as: 

A L (g)=L- 1 f(gLi) , (4) 

where f(gL^) is a dimensionless scaling function and v is 
the critical index of the correlation length. The unknown 
scaling function / is determined by the best collapse of 
the finite size ED data. The fitting gives v w 0.75 (cf. 

Fig. p. 



B. String Order Parameter 

In the absence of the conventional Landau order pa- 
rameter, the different gapped phases in the ladder can 
only be distinguished by a non-local string order param- 
eter (SOP), which is non-zero in one phase and vanishes 
in the other. The SOP was first introduced in 1989 by 
den Nijs and Rommelse^ who related it with the gapped 
or Haldane phase in the spin-1 Heisenberg chain. In that 
system there is no local symmetry breaking and therefore 
no local order parameter can be identified. However, the 
SOP detects the hidden ® Z>2 symmetry breaking^ in 
the Haldane phase. The SOP notion was later general- 
ized for two-leg spin ladders asJ^~— 
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FIG. 4: Two- leg ladder; critical line J± C (S) where the gap of 
the staggered phase vanishes. Adapted from Ref. [[1, original 
data from Ref. @. 




FIG. 5: (Color online) The collapsed scaling function (left panel) and gap (right panel) for the staggered two-leg ladder with 
S = 0.5, calculated for N = 8 through 28 along with the fit using optimized parameters v ~ 0.755 and J± c w 1.27. 



m— 1 

°odd/eve n = - J"* (^ung/dmgM eX P J2 5 rung/diag(0 }S£ing/dtag("0 ) . ( 5 ) 

J=n+1 



where S r * ung ( n ) = Sf (n) + Sf (n) and Sf iag (n) = Sf(n) + 
5|(n + 1). The even and odd SOP have been shown to 
be mutually exclusive in several numerical and analyti- 
cal studies . 11 ' 23 Previous numerical studie o 11 ' 19 have as- 
sumed a staggered dimerization, while here the SOP were 
calculated for both configurations by the ED method. In 
agreement with the previous studies, we find that the 
SOPs are featureless for the columnar phase: in the whole 
phase space domain (J±,5) the even string order param- 
eter is non-zero, while the odd string order parameter 
vanishes. 



We find for the staggered phase that the even/odd SOP 
is non-zero only in the rung-dimer/leg-dimer phases, re- 
spectively, and vanishes at the critical point. The calcu- 
lated SOPs are shown in Fig. [6] The SOP vanishes at the 
critical point as OJ oc g 2/3 where j3 is the critical index of 
the order parameter ! 20 ! 26 Then the SOP in the system of 
length L scales as 

OXg) = L-%f(gL$) , (6) 

where the scaling function / is determined as for the 
gap by the requirement to collapse the numerical data 
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FIG. 6: (Color online) The odd and even string order param- 
eters of the staggered two-leg ladder with 8 — 0.5. The ED 
data for ladders with number of spins N = 12 to 28 are shown, 
along with the critical scaling curves. 



TABLE I: Critical scaling indices of the two-leg ladder. 
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even/odd 
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Jxc/J 


0.25 


even 


0.05294 


0.758 


0.775 


0.25 


odd 


0.285 
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even 


0.0251 


0.726 


1.32 


0.5 


odd 


0.215 


0.698 


1.21 


0.75 


even 


0.0687 


0.764 


1.62 


0.75 


odd 


0.235 


0.648 


1.64 


1.0 


even 


0.0124 


0.7306 


2.08 


1.0 


odd 


0.227 


0.7013 


1.98 



for the systems of different sizes onto the same curve. 
From the analysis of the ED data for S = 0.5 we obtain 
a critical index /3 dd ~ 0.25 for the odd string order and 
/3evcn ~ 0.025 for the even string order. The collapse of 
the scaling functions are shown for 6 = 0.5 in Fig. [7] The 
same critical scaling analysis was performed for several 
values of 5, and these results are tabulated in Table [H 



III. THREE-LEG LADDER 



Similar to the two-legged case the dimerized three-leg 
ladder occurs in two distinct configurations, the stag- 
gered and columnar patterns (cf. Fig. [5]). The Hamil- 
tonian of the intrinsically dimerized three-leg ladder is 
given by Eq. ([T]) with m = 3 and like the two-legged 
ladder the staggered/columnar patterns are given by the 
coupling J a (n) defined by Eqs. ©/CI]) respectively. 



A. Ground State Energies and Gaps 

The ground state energies were calculated by the ED 
method for the three-legged ladder with N = 24 total 
number of spins. The columnar pattern was found to 
have the consistently lower energies throughout the pa- 
rameter space, in agreement with the previous results 
from the bond mean- field theory^ The ED results for 
the energy gaps, given in Fig.[S]for N = 24, show that the 
columnar configuration is always gapped, and the stag- 
gered pattern approaches gaplessness at a critical value 
of rung coupling J± c . Similar to the two-legged case, the 
appearance of criticality in the three-legged ladder de- 
pends on the dimerization pattern. The phase diagram 
of the staggered ladder is given in Fig. [KjJ The ground 
states of these ladders have been found to be well de- 
scribed by short-range valence bond solid (VBS) states. 
A particular VBS state is denoted according to the num- 
ber of valence bonds formed with contiguous states; i.e., 
(m, n)-VBS with m + n = 35. A transition from a (1,2)- 
(VBS) phase to a (2,1)- VBS phase is understood to occur 
on crossing the phase boundary^ The finite size data 
for the energy gap was used to extract the critical in- 
dices near the transition point. This was done using a 
critical scaling analysis similar to that used to study the 
two-legged case, with the critical scaling relation given by 
Eq. ((4]). Data for laddders with N = 12 to 30 were used. 
The scaling exponents on each side of the critical point, 
denoted as Region (1,2) (g < 0) and Region (2,1) (g > 0) 
were found to be different, with separate values ^(1,2) 
and ^(2,1) ■ Collapse of the scaling function for 5 = 0.5 
(cf. Fig. |I2J yields i/ (1>2 ) » 0.863, v {2 ,i) ~ 0.385 and 
J± c w 1.08. For (5 = 1, the fitting gives v^ 2 ) » 0.789, 
^(2,1) ~ 0.363 and J± c ~ 1.65. As opposed to the known 
behavior for large values of rung coupling J± of the gap 
of the two-legged ladder, which increases as A oc J±, the 
gap of the three-legged ladder quickly saturates in both 
dimerization patterns (cf. Fig. [TTT) . This fact limits the 
width of the critical region where the finite-size scaling 
analysis is valid. 



B. String Order Parameter 

The generalized SOP for three-leg ladders is given by 
Eq.© where 5 r z ung (n) = Sf(n) + 5|(n) + S§(n) and 
S z diag (n) ee 5f(n) + Sl(n + 1) + Sf (n + 2). The SOP 
were calculated for N = 24 for both the staggered and 
columnar configurations. Results for smaller systems (N 
= 12-18) give unreliable results due to the very small 
number of rungs in the sum. For this reason, a critical 
scaling analysis can not be performed for an estimate of 
j3 until much larger system sizes (N = 30-36) are reached. 
Similar to the two-legged case, the columnar pattern was 
found to possess even string order throughout the param- 
eter space, with no critical features. The staggered phase 
of the ladder possesses even string order associated with 
the (1,2)-VBS phase for < S < 5 C , which vanishes upon 
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FIG. 7: (Color online) The lefthand (righthand) panel shows the fitted scaling function for the odd (even) string order parameter 
of the staggered two-leg ladder with S = 0.5. The data collapse for the optimized scaling parameters v — 0.7, J± c = 1.21, and 
/3 odd = 0.25 (u = 0.74, Jxc = 1.32, and /3 cven = 0.025). 
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FIG. 8: Dimerized three-leg ladder. Bold/thin/dashed lines 
represent the stronger/ weaker chain coupling J(l ± S) and 
rung coupling J±, respectively. Dimerization patterns: (a) - 
staggered; (b)- columnar. 
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FIG. 9: (Color online) Energy gaps for the staggered and 
columnar three-leg ladder dimerization patterns calculated by 
the exact diagonalization for N = 24. 



entering the (2,1)-VBS phase, accompanied by emerging 
odd string order. Representative results for the SOP in 
the staggered ladder are given in Fig. I14I However, we 
note that a residual odd string order persists in the (1,2)- 
VBS region of the parameter space. We believe this effect 
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FIG. 10: Three- leg ladder; critical line J± C (S) where the gap 
of the staggered phase vanishes. Adapted from Ref. Qj| orig- 
inal data from Ref. [ill. 
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FIG. 11: (Color online) Three-leg ladder; gap saturation for 
strong rung couplings J±. 
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FIG. 12: (Color online) The collapsed scaling functions for the staggered three-leg ladder with 5 = 0.5 in Region (1,2) (left 
panel) and Region (2,1) (right panel), calculated for N = 12 through 30 with optimized parameters f(i,2) ~ 0.863, ^(2,1) ~ 0.385 
and J ±c « 1.08. 
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FIG. 13: (Color online) The ED data for the ground state 
energy gap with 8 — 0.5 for number of spins TV = 12 to TV = 
30, along with the fit using optimized parameters fn.g) = 
0.863 in Region (1,2), f(2,i) = 0.385 in Region (2,1) and J± c = 
1.08. 



is due to the very small size of the system, and should 
disappear if larger system sizes are assessed. 



IV. CONCLUSIONS 

We analyze the dimerized antiferromagnetic spin-i 
two- and three-leg ladders. In agreement with the earlier 
mean-field predictions,— the present ED results demon- 
strate that the columnar phases have lower ground state 
energies than the staggered phases. Thus, if no mecha- 
nism for locking into a particular dimerization pattern is 
provided, a dimerized ladder would order in the colum- 
nar phase. In this paper we analyze the two possible 



dimerization patterns on the same footing, as intrinsic 
given properties of the ladders without providing specific 
mechanisms responsible for a particular dimerization. 

The existence or absence of the quantum critical point 
is dependent on the dimerization pattern. The gapped 
phases, which cannot be distinguished by the conven- 
tional Landau long-range order parameter, have been 
characterized using the non-local SOPs, which are non- 
zero in one phase and vanish in the other. Here, for the 
first time, we systematically identify and numerically cal- 
culate all SOPs for the both critical (staggered) and non- 
critical (columnar) phases in the two-dimensional space 
of model's couplings (S,J±/J). In addition we analyze 
the critical parameters near quantum criticality. 

For both 2- and 3-leg ladders the columnar pattern was 
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O z 




FIG. 14: (Color online) The odd and even string order param- 
eters of the staggered two- leg ladder with J±/J = 0.25 and 
J±/J = 0.75 as functions of dimerization S. The ED data for 
ladders with N = 24 total number of spins is shown. 
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found to be always gapped and to possess the non-zero 
even SOP which is essentially featureless. The finite-size 
ED data suggest that the odd SOP vanishes in the limit 
N — > oo in the columnar phase for both types of ladders 
considered. 

The finite size data for the energy gap were used to 
extract the critical indices in the vicinity of the line of 
quantum phase transition in the staggered case. The 
scaling behavior of the three-leg ladder was found to be 
different on each side of the critical point. These results 
give for 6 = 0.5 values of ^(1,2) ~ 0.863 and f(a,i) ~ 0.385. 

From the critical scaling analysis of the ED data for 
the ladder with two legs we find the critical indices of the 
correlation length v a 0.755, of the odd order parameter 
/5 dd ~ 0.2 and of the even order parameter of /3 eV en ~ 
0.02. The rest of the critical indices can be determined 
by the scaling relations. 

Most likely, the smallness of /3 C vcn and large difference 
between /3 evon and /3 dd are due to stronger finite-size ar- 
tifacts in 0^ vcn . The non- monotonous behavior of 0* ven , 
clearly seen in Fig. [51 also makes the finite-size scaling 
analysis less reliable. 

According to the predictions of Refs. the crit- 

ical behavior of the two-leg ladder along the critical 
line J_L C (^) is the same as near the integrable point 
(J±/J = 2, 8 = l). 27 Since along the line S — 1 the ladder 



([T]) reduces to a dimerized Heisenberg chain, the above 
assumption implies a critical behavior with v = 2/3 and 
/3 = 1/12 modified by multiplicative logarithmic correc- 
tions due to the marginally irrelevant term in the effective 
Hamiltonian. Our values v sa 0.75 and j3 « 0.2 (including 
the case 6 = 1) are comparable with these predictions, 
especially taking into account: (i) natural limitations on 
the precision for the critical indices due to the small sys- 
tem sizes which can be reached by ED; (ii) fitting the 
critical indices without logarithmic corrections exagger- 
ates their values (see, e.g., [28.29]). The critical behavior 
of the ladders in the vicinity (along) of the critical lines 
clearly warrants a detailed analytical study. We plan to 
address this issue in our future work. 
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